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Abstract 



We derive and solve kinetic equations for leptogenesis within the Closed Time 
^ ■ Path (CTP) formalism. It is particularly emphasised how the procedure of real 

intermediate state subtraction familiar from the Boltzmann approach is realised 
within the CTP framework; and we show how in time-independent situations no 
lepton asymmetry emerges, in accordance with the CPT theorem. The CTP 
approach provides new quantum statistical corrections from evaluating the loop 
integrals. These lead to an enhancement of the asymmetry originating from the 
Bose statistics of the Higgs particles. To quantify this effect, we define and evaluate 
an effective CP-violating parameter. We also solve the kinetic equations and show 
explicitly that the new quantum statistical corrections can be neglected in the 
strong washout regime, while, depending on initial conditions, they can be sizable 
for weak washout. 



1 Introduction 



In the light of the present knowledge of elementary particle physics phenomenology, lep- 
togenesis appears as one of the most plausible explanations for the dynamical emergence 
of the baryon asymmetry of the Universe [H [2]. Leptogenesis in turn is a variant of 
baryogenesis from out-of-equilibrium decays of heavy particles during the expansion and 
cooling of the early Universe [31 H] . The deviation from equilibrium provides the break- 
ing of time reversal invariance, that is a necessary condition for baryogenesis j5] as a 
consequence of the CPT theorem. 

The most common approach for making quantitative predictions is to formulate and 
solve a network of Boltzmann equations, that describes the evolution of the number 
densities of the particle species relevant to baryogenesis. In simple scenarios for leptoge- 
nesis, a setup that leads to meaningful results is to compute the distribution function of a 
heavy right-handed Majorana neutrino Ni and of the charge density of a lepton doublet 
£, while one may assume that the remaining species are in thermal equilibrium. More 
realistic models typically take account of the presence of the several generations. The 
Boltzmann equations are usually written in a form, where on the left hand side, there 
is a kinetic term that traces the particle number distributions, while on the right hand 
side, there is a collision term. The collision term accounts for elastic and inelastic (par- 
ticle number violating) scatterings. It arises from weighting reaction rates, computed 
in vacuum scattering theory (in-out formalism), with the distribution functions of the 
species involved. 

When setting up the Boltzmann equations, one encounters the following issue: Let 
rVi-^, the rate for the decay process N — > £<fi, be proportional to 1 + s, where e 
parametrises the CP asymmetry of the decay. Then, by complex conjugation of the 
involved coupling constants, it follows that T Ni ^ f i ( j } * oc 1 — s. Application of the CPT 
theorem then implies that the inverse decay rates Te^Ni oc 1 — s and r^*^^ oc 1 + s. If 
these were the only contributions that are accounted for in the collision term, one might 
conclude that the two-by-two rate £<p — > £<ft* is proportional to 1 — 2s, while the rate of 
£<p* — > £(j) is proportional to 1 + 2s. An asymmetry even in equilibrium would result, in 
apparent contradiction with CPT invariance. But of course, there is no contradiction, 
since the heavy neutrinos are unstable and therefore do not belong to the Hilbert space 
of asymptotic states of a unitary S'-matrix. The total lepton number violating two- 
by-two rates are therefore not given by the above naive multiplication of the rates for 
decay and inverse decay. Yet, it is crucial that one tracks the neutrino densities in the 
Boltzmann equations. The issue is fixed when realising that the resonant parts of the 
two-by-two scatterings, where the intermediate neutrino is on-shell, contribute to the 
Boltzmann equations at the same order in coupling constants as the decay and inverse 
decay processes. Technically, one adds the full two-by-two rates and subtracts from these 
the on-shell contribution, since this is already accounted for in the decays and inverse 
decays [I]. This procedure is known as the subtraction of real intermediate states (RIS). 

One expects that the RIS subtraction is naturally realised within a first principle 
approach, that is less heuristic than the Boltzmann ansatz. This first principle approach 
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is given by applying the Schwinger-Keldysh Closed Time Path (CTP) formalism [61 [7] 
to field theory [8j. In this paper, we derive and solve the Kadanoff-Baym equations 
for leptogenesis, that result from the CTP formalism. We pay particular attention to 
demonstrating that no lepton asymmetry is generated in equilibrium, or, more generally 
in situations that are time-reversal invariant. For models with fermionic leptons this 
result has been reported earlier in works based on the CTP formalism [9KT011II], though 
no detailed derivation has been given. (The main focus of Refs. [TUJ [11] is the study 
of memory effects.) For the case of scalar leptons, the vanishing of the asymmetry in 
equilibrium has been demonstrated recently in Refs. [I2J H3] • 

In addition to offering the framework for a first-principle derivation of equations 
that describe leptogenesis, the CTP formalism also provides finite density corrections 
that do not straightforwardly follow in the Boltzmann approach. These corrections take 
the form of an integral over terms linear in the Higgs boson and lepton distribution 
functions. Historically, it has been proposed to compute the CP asymmetry from the 
discontinuity of the correlation functions in the Matsubara (imaginary time) formalism. 
In the CTP (real time) approach, this corresponds to the imaginary part of the time- 
ordered correlation functions [HI [15]. It is characteristic of these Green functions, that 
they feature not only terms that are linear in the lepton and Higgs boson distribution 
functions, but also a product term of these distributions. In contrast, it appears to be the 
consensus among the more recent papers [TUl [TT| [T2 l ITS] that the CP asymmetry does not 
receive corrections from a product term of Higgs boson and lepton distribution functions. 
This has been shown in detail for the model with scalar leptons in Refs. [T2| [T3] . Here 
we find the same property in the model with fermionic leptons: The CP asymmetry 
results from terms that can be written as the difference of Wightman type correlators 
(no explicit time-ordering of the field operators). 

After deriving the kinetic equation for the lepton asymmetry in the CTP formalism 
we solve it numerically in various situations. We find that the effect of the new quan- 
tum statistical corrections is small when the bulk of the lepton asymmetry is produced 
at temperatures T <C Mi, where Mi is the mass of the neutrino that decays out of 
equilibrium, because the distribution functions are strongly Boltzmann suppressed for 
momenta much large than T. This situation corresponds to the strong-washout scenario, 
which is of particular interest since the resulting asymmetry is independent of both, the 
initial abundance of heavy neutrinos and a possible primordial lepton asymmetry. On 
the other hand, when the bulk of the asymmetry is produced for T ~ Mi or larger, 
which corresponds to weak washout, the impact is more sizable. For a vanishing initial 
lepton asymmetry and thermal initial conditions of N\, we find an enhancement of the 
asymmetry, as apparently the Bose enhancement from the Higgs particle distribution 
function dominates over the Fermi suppression from the lepton distribution. This qual- 
itatively confirms what has been reported for scalar leptons [T2"| [13]. In this context, we 
note that the effects of finite density have been considered when the quantum statistical 
factors are accounted for in the initial and final states [TH1 EE] , but not in the loop cor- 
rections, in which case it has also been found that the corrections may be sizable in the 
weak-washout regime. However, the quantum statistical factors in the loop corrections 
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appear within the CTP approach on the same footing and with the same size as the 
external state factors, which is why they cannot be neglected consistently. 

The remainder of this paper is organised as follows: In Section [21 we specify the model, 
introduce the free CTP propagators and provide the one-loop self energies. These ingre- 
dients are used in Section [3] to derive the Kadanoff-Baym equations that correspond to 
the tree-level Boltzmann equations with quantum statistical factors. The wave-function 
correction that leads to a CP asymmetry arises from a two- loop diagram in the CTP 
formalism, which we discuss in Section HI We explicitly show that no asymmetry is pro- 
duced in equilibrium and we relate this fact to the usual procedure of RIS subtraction. 
Similarly, we discuss the vertex correction in Section In Section El we numerically 
evaluate the effective CP-violating parameter and solve the Boltzmann equations for 
different strengths of washout, which demonstrates that finite density corrections are 
important for weak washout, but may be neglected in the strong- washout domain. We 
conclude in Section [71 

2 CTP Approach to Leptogenesis 

In the simplest scenario of leptogenesis, there are two right handed neutrinos iVj (i — 1, 2) 
with lepton- number violating Majorana masses Mj, one scalar SU(2) Higgs doublet <fi 
and an active left-handed lepton doublet £. The Lagrangian of the model is 

C = - Mi)ip Ni + i>oW>t + {d^){d^) - Y*^P R i/j m - Y^mP^t , (1) 

where a summation over % = 1,2 is understood and -Pl/r = (1 -F 7s)/2- Besides, it is 
implicitly assumed that the SU(2) indices are contracted in a gauge-invariant way, that 
is (fiipi = —(pA^AB^tB and -ip^ = (^)a^a_b(0^)_b, where €ab denotes the antisymmetric 
2x2 matrix with eu = 1. Upon the four- component spinor ipm, we impose the Majorana 
constraint 

1>h = Cffii = iim , (2) 
with C the charge conjugation matrix. 

2.1 Propagators 

For the CTP formalism we follow closely the notational conventions of Refs. [T81 [TOl [20] 



and define the CTP propagators for the leptons by 

iS eap( u ' v ) = iS f a pi u > v ) = -(<M<#fo(u)) , (3a) 

i5 S( u > v ) = iS i a p( u > v ) = (i>ia{u)$tp{v)) , (3b) 

iS £J( u > v ) = iS Lp(u, v) = (T{il)f, a {u)$w{v))) , (3c) 

iS eap( u i v ) = iS iap(u, v) = . (3d) 
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We perform a perturbation and gradient expansion of the full Kadanoff-Baym equations 
that uses the tree-level propagators and vertices as elementary building blocks. The tree- 
level propagators are therefore given by the vacuum propagators plus a finite-density 
contribution proportional to the spectral function of the free theory. When there are 
no correlations between modes of different momentum and spin, we can characterize 
the state by particle number distribution functions. In case there is also no preferred 
direction for the spin, the usual averaging procedure applies. These assumptions lead to 
the explicit expressions 

iS<(p) = -2n5{p 2 )P^P R [0(po)fi(p) ~ 0(-po)(l - Fi(-P))] , (4a) 
iS>(p) = -2n5(p 2 )P^P R [-0(p«,)(l - fi(p)) + H~Po)M-p)] , (4b) 

= Pl ^T^ Pr " 27r5 ^ 2 ) PL ^ [#(Po)ft(P) + V(-Po)M-v)] , (4c) 

LSf (p) = -P L ^-P R - 2n5{p 2 )P^P R [0(po)fi(p) + #(-Po)M-p)] , (4d) 
p z — is 

for the lepton propagators, where fe(k) denotes the distribution function of leptons 
and /g(k) of anti-leptons. The leptons i occur within an SU(2) doublet, but since this 
symmetry is unbroken at the time of leptogenesis, we assume that the lepton densities 
are SU(2)-symmetric and hence diagonal. In particular, the lepton propagators defined 
here are the diagonal components of a diagonal propagator with additional SU(2) indices: 

S™j 2 \ u ,v) = 5 AB S e (u,v) A, B = 1,2. (5) 
Similarly, for the Majorana-fermion propagators we find 

iSUp) = ~^5(p 2 - M 2 ){j> + Mi) [0(po)/jw(p) - 0(-Po)(l " fm(-p))] , (6a) 
iS^P) = ~^(p 2 ~ M 2 )(i> + Mi) [-7?(po)(l - fm(p)) + H-Po)fm{-p)} , (6b) 

iS m(P) = ^tff- " ^(p 2 - M 2 ){i> + Mi) [#(po)fm(p) + #{-po)fm{-p)] , 

p z — Mj + IS 

(6c) 

= - J-te-ir - ~ + M ^ VMM + V(-P0)fm(-P)] ■ 

p Ivl^ Lb 

(6d) 

As a consequence of the Majorana condition (j2J), the distribution functions for neutrinos 
and anti-neutrinos are identical, and the propagators inherit the property 

S Ni (u,v)=CS t m (v,u)C\ (7) 
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where the transposition acts here on both the CTP and the Dirac indices. Finally the 
scalar propagators take the form 

iAj(p) = 2n5(p 2 ) [0(po)/*(p) + ^(-Po)(l + /*(-p))] , (8a) 

iA>(p) = 27r5(p 2 ) [tf(p )(l + U(p)) + ^(-Po)^(-P)] , (8b) 

iAj(p) = + 2vr5(p 2 ) [<?(po)/«(p) + ^(-Po)/0(-p)] , (8c) 

= + 27r ^ [*(p«)Mp) + *(-Po)/*(-P)] • ( 8d ) 

As for the leptons, it is understood that 

A^gW) = WV^) A B = 1,2. (9) 

In the following, an occasional multiplicity factor g w — 2 arises from performing SU(2) 
traces such as 

tr su{2) [A^e+Sf (2) e] = A^ tr su(2) [e+e] = g w A^S e . (10) 

In thermal equilibrium the distribution functions assume the usual Fermi-Dirac and 
Bose-Einstein forms, respectively: 

MP) = ft(p) = fP(p) 



0\P\ + 1 



e 



A»(p) = /j%(p) = a r^- 2 , (11) 

e P\/P 2 +Mf _|_ 1 

/*(p) = /*(p) = /; q (p) = • 

Furthermore the equilibrium Green functions satisfy the Kubo-Martin-Schwinger (KMS) 
relation, 

iS>(p) = -e^iS<(p), iS> l (p) = -^iS< i (p), iA>(p)=e*»iAJ(p), (12) 
an extremely useful property that we employ exhaustively throughout this paper. 

2.2 One-loop self-energies 

The one- loop neutrino self-energy in the CTP formalism is given by (cf. the vacuum 
expression in Ref. |21j ) 

vffivW =g w j 0^0^^ w ( k ~ k '- k ") (1 3 ) 

x [YiYjPdS? {k')P K iAf(k") + Y*Y 3 C [P L iS b e a (-k')P R ] T CHA%(-k") 



where a, b — ± denote the CTP indices, see Eqs. fl3aj) - fl3dl) . Our definition for the self- 
energies is such that — iE corresponds to the sum of all one-particle irreducible diagrams. 
The one-loop lepton self-energy reads 

in " {k) = |Fi|2 / W^WT^^^ ~ k '~ ^^OaiA^-fc") , (14) 

where we sum over the neutrino flavours i. In thermal equilibrium these self-energies 
inherit the KMS property 

Vm/p) = -e^ • (15) 

from the constituent CTP Green functions. At the one-loop order it is easy to verify this 
explicitly by inserting the tree-level propagators (BJ, ([6]), (IE]) into the self-energies ( TT3l) . 
( JT4|) . One of the main objectives of this paper is to show explicitly how KMS also remains 
valid when considering the next-to-leading order contributions to S^v^. It turns out that 
this encompasses the subtraction of RIS as it is performed in the conventional approach 
to baryo- and leptogenesis. 

3 Tree-Level Contributions to the Kinetic Equations 

The theoretical description of leptogenesis is concerned with the time evolution of non- 
equilibrium densities of active leptons and of right-handed neutrinos. This is described 
by the Kadanoff-Baym equations (see e.g. [1HJ [X9l |20] ) . We perform the standard Wigner 
transformation and gradient expansion. Since we are interested in finite-density effects 
from loops in the CTP formalism, we restrict ourselves to the leading term in the gradient 
expansion. We regard this as a self-consistent leading-order approximation to the full 
Kadanoff-Baym equations, relative to which corrections from higher gradients and loops 
may eventually be considered. Thus, the kinetic equations take the form 

\i^(k)P L iS<(k) - i$f(k)P L iS>(k)] , (16a) 

[i9i«(A0LS&(A0 - i?$«(*)LS&(A0] , (16b) 

where we recognize the familiar expression for the collision term in the CTP formalism 
on the right-hand side. Here and in the following we assume spatial isotropy, which 
means that there is no angular dependence in the distribution functions and we may 
identify ft^mQ*) = ft,<i>,m(M). In particular, we repeatedly use fe^ t m(-^) = fi^mfe) 
to simplify the arguments of distribution functions. When performing the Wigner trans- 
formation and gradient expansion in a slowly varying background, all time dependence 
in the propagators (j3J, (jSJ), (JE} is isolated in the distribution functions fe^m, where we 
suppress an explicit time argument. The equation for the lepton density can be simpli- 
fied when taking the Dirac trace of Eq. (|16ap . multiplying by minus one and integrating 
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over fc . Inserting Eq. @, we obtain 

j t (f £ (k) - /i(k)) = C<(k) = | ^tr l$f>{k)PvSf{k) - iVf(k)P L iS>(k)] , (17) 

which defines C^(k). In a similar way we obtain the equation for the time dependence of 
the Majorana neutrino distribution function 

j t fm(V = C m (k) = \j ^sign(fc )tr [^(k^k) - i^M^k)] . (18) 
3.1 Tree- level collision terms 

The tree-level collision term for the leptons can be straightforwardly evaluated by insert- 
ing the tree-level propagators into Eq. (ITTj) and the expression (j!4p for the self-energy. 
Since we are interested only in the lepton number density 

^ = Sw? m ' (19) 

we integrate Cg(k) over k and find 

; w = -in* [ ( 2 °) 



(2tt) 3 " y ' 1 11 J (27r) 3 2|k 2 | (27r) 3 2v/k' 2 + M? (2vr) 3 2|k" 

x (2n) i 6\k l - k - k") 2k' ■ k |(1 - f Ni (k')) x [/<(k)/,(k") - /*(k)^(k")] 
- /^(k') x [(1 - fl(k))(l + /^(k")) - (1 - £(k))(l + / (k"))] 



Due to the fast violation of Higgs number, we can assume that — f<j> = f^ q . For the 
leptons we may assume fe — = —(fe — fl q ), as justified by fast pair-annihilating and 
-creating interactions among the leptons. Expanding to linear order in the deviations 
from equilibrium, the lepton collision term simplifies to 

r d 3 k _ 2 r d 3 k d 3 k' d 3 k" 

J (2vr) 3 H J " il J (27r) 3 2|k| ( 27 r) 3 2 x /k^TMf (27r) 3 2|k"| 1 ] 



x 



(27r) 3 2|k| (27r) 3 2v/k' 2 + A^(27r) 3 2|k" 
{2^)H\k> - k - k") 2k' ■ k [f;\k") + f m (k')] x [f e (k) - f £ (k 



This term describes the washout of a lepton asymmetry through inverse decays. In a 
similar way we obtain the tree-level collision terms for the neutrinos by inserting the 
tree-level propagators into Eq. (ITS]) and the expression (fill) for the self-energy. For 
the neutrinos, we shall need the phase-space distributions, so we need the unintegrated 



7 



collision term, for which we find 

c - (k) 4^' 2 ^hm I (27r w * - fc/ - ° (22) 

xtr [>'(# + Mi) {-/j«(k)[l - /,(k')][l + U(k")} + [1 - f Ni (k)]Mk')U(k")} 

+ Mi) {-/j«(k)[l - fe(k')][l + U(k")} + [1 - /^(k)]/,(k')^(k")} 

It is useful to check that when setting the Higgs and lepton distribution functions to 
zero, this reproduces the negative total decay rate of Ni, 

--fam-jg^fM. (23) 

Next, introducing 

Sfm = fm - fZ (24) 

and expanding to linear order around equilibrium distributions the Majorana collision 
term simplifies to 

1 f d 3 k' d?k" 



C m (k) = - gJYA 2 — - = / - " , - ^ , (27r) 4 (5 4 (fc — k' — k") (25) 

1 ^ ' ' 2v/k 2 + Mf J (2vr) 3 2|k'|(27r) 3 2|k"| l ; 1 1 K ' 

x2k.k'Sf Nl (k)[l-fP(k') + f;\k")]. 

Note that SJni is not necessarily small in some realistic scenarios for leptogenesis, e.g. 
for weak washout. 

It is easy to see from Eqs ( 12 ip . ( 12 5 j) that in thermal equilibrium the collision terms 
vanish. Alternatively, one may either by directly use the KMS relations ( JT2lfT5|) in 
Eqs. (j!7|18p or explicitly insert the equilibrium distributions into Eqs. ( 120|22|) . This 
means that the equilibrium state is a static (time-independent) solution to the Kadanoff- 
Baym equations. 



3.2 Expansion of the Universe 

Before proceeding we comment on how the effects of particle dilution due to the ex- 
pansion of the Universe can be systematically incorporated into our equations. So far 
we have formulated the Kadanoff-Baym equation on the Minkowski background. Ne- 
glecting Planck-scale suppressed corrections the only effect of replacing the Minkowski 
background by the Friedmann-Robertson- Walker one is that the particle modes are red- 
shifted. 

This can be conveniently implemented by adopting conformal coordinates, which are 
defined by the metric 

SV = a 2 (v) diag(l, -1, -1, -1) , 
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where t] = x° denotes the conformal time and a(rj) the scale factor. Then the La- 
grangian ([1]) describes the dynamics in the expanding background, provided that we 
have redefined the fields such that their kinetic terms are canonically normalized, that 
we replace Mj — > aMi and that we assume that the Higgs field has no coupling to the 
scalar curvature R. Consequently all expressions given earlier in this Section and in Sec- 
tions H] and [5] below remain valid in the expanding Universe if momenta k are understood 
as comoving momenta k com , temperature T (and (3 = 1/T) as comoving temperature T com 
(and /3 com ), time t as conformal time rj, /(k com ) and n as comoving phase-space distribu- 
tions and particle number densities, respectively, and if finally masses M, are replaced 
by aMi, wherever they appear. 

For the scale factor in the radiation-dominated Universe we use the relation 

a{v) = or 7 ? • (26) 

Using the Friedmann equation and the free energy, we can relate this to the temperature 
as 

T= i_ i /^fjs y /4 , (27) 



a(r]) a(rj) \ 2 V5'* 71 ' 3 . 

where denotes the number of relativistic degrees of freedom and mpi = 1/ \/Gn = 
1.22 • 10 19 GeV the Planck mass. The variable z = M\jT that will be used in Section [6] 
is therefore related to rj by a constant proportionality factor. 

Kinetic equations are often formulated in terms of physical momenta and physical 
time. In order to make contact, we can replace in our equations the comoving momenta 
by physical momenta /c p h = k com /a and the comoving temperature by the physical tem- 
perature T p h = T com /a. Doing so and dividing by a factor of a for convenience, the 
right-hand sides (i.e. the collision terms), simply take the same form as given earlier in 
this Section, where now all momenta are to be understood as physical. The left-hand 
sides change according to 

* /(k com ) = Jr/(V) + ( J-f(k ph )) %J- (28) 



a(r)) dr) dt \c?|k ph | / dt 

= |/(k ph )-i/|k ph |^ / (k ph ), 

where / may stand for Jn or ft— ft. This formula applies to the treatment of the densities 
of right-handed neutrinos, that do not necessarily maintain kinetic equilibrium during 
leptogenesis, implying that one cannot substitute for f^i a Fermi-Dirac distribution 
with a (pseudo-) chemical potential. For the lepton charge, it is suitable to consider the 
integrated version 

d% h ^L /(kcom) = d + 3Hn ^ (2g) 



a(ri) J (2tt) 3 dr] J K dt 

of the above equation. This is just the familiar cosmological dilution law, where n p h 
denotes a physical number density and H the Hubble parameter with respect to time t. 
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However, in Section [6] we will solve the kinetic equations directly in conformal coordi- 
nates. 



3.3 Boltzmann equations 

We now write down explicitly the Boltzmann equations that describe the time evolution 
of the distribution function f^i and of the lepton asymmetry in conformal coordinates: 

■jj-fmiKaw) = D(k com ) , (30) 

^(n e -n e ) = W + S. (31) 

Here D(k com ) denotes the contributions ( l25i) from decays and inverse decays of N±, 
while W denotes the washout term |UJ). S = S wi + S v is the source term for the lepton 
asymmetry, which we decompose into a wave-function contribution S and a vertex 
contribution S v . Note that here and below, unless indicated otherwise, rig refers to the 
comoving particle number density. 

Performing the substitution of variables pertinent to conformal coordinates as dis- 
cussed above in Eq. (1221) . we obtain: 

D(k com ) =2 9w \Y^-^==L=== / , "t- l J" c Z , (32) 



2 v /k c 2 om + a2(r / )M 1 2 i (27r)32|k' com | (27r)»2|k» 
x (27r) 5 (k cora — k com — k com )k com ■ k cora 

- /7vi(k com ) (1 - fa (k com ) + (k com )) + f e (k com ) (k com ) 

In the evaluation of the tree-level collision term we may set fg and f^ equal to the equi- 
librium distributions given in Eq. (ITT]) (but expressed in terms of comoving momentum 
and temperature), which allows us to evaluate the phase-space integrals. It proves useful 
to define the quantity 

^ = ^j ( 27r )32| P | (27r ) 3 32| q | ^ 4 ^ k -*> - g)*>» (* - /Hp) + ;; q (q)) , (33) 

which we will calculate analytically in Section |6] In terms of this we find for the decay 
and inverse decay term 

2 k com r c i 

-D(kcom) = — 2|Yi| — ^ N {k com ) /Ari(k com ) /jvi(kcom) ) (34) 

where fc com0 = \fk 2 cora + a\j] 2 Ml and 

/^ q i( k com) = efccomo/rcom + I ■ ( 35 ) 
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To calculate the washout term (I2ip we approximate f^i by the equilibrium distribution 
and assume that the leptons are in kinetic equilibrium so that fi can be parameterized 
in terms of a chemical potential. Expanding to linear order in the chemical potential, 
the latter assumption yields the relation 

f e (k) - f e (k) = (n, - n,) x £ 2 . (36) 

(ePi k i + 1) 



We then obtain for the wash-out term 



d 3 k „ „ \ 31yd 2 a 2 M, 2 , _ , f , e x 

^ = / Ce{k) = ~^t^ {ni " nd I dx W+W ln 



/ a 2 M, 2 \ 



e 4 ^m^ ' " + 1 

a 2 M 2 



(37) 



The substitution f Ni — > which was made to take analytical integrations further, 
is strictly speaking only a correct approximation when Sf^i As in scenarios 

of strong washout or in the transition region from weak to strong washout most of 
the asymmetry is eventually produced when Aq is close to equilibrium, the error from 
assuming an equilibrium abundance of N% in the washout term at all times should be 
small in the final result. (In fact, when assuming Maxwell statistics [2J [T5], [22] the absence 
of Pauli blocking factors in Eq. ( 1201) implies an even more drastic approximation, since 
then f Ni (k') -> in Eq. (I2T]>.) 

As is well known, the source term S for the lepton asymmetry in Eq. (|3T|) only arises 
in fourth order in the couplings Yj. Deriving this source term including all quantum 
statistical factors is the subject of the following two sections. 



4 Self-Energy Contribution to the CP Asymmetry 

A key ingredient to baryogenesis calculations in the conventional framework is the sub- 
traction of RIS. If this were omitted, a baryon asymmetry would result in thermal equi- 
librium in contradiction with the CPT theorem and the unitarity of scattering matrices. 
In this Section we show how the RIS subtraction arises within the CTP framework for 
leptogenesis and derive the CP asymmetry from the wave-function correction in the 
presence of finite number densities. 

4.1 Wave-Function Correction to Cg 

The aim is the calculation of the wave-function correction to the lepton collision term 

(ED, 

Cf(k) = [ ^tr fe7 f> (fc)P L iS<(A;) - ^(kWSUk)} , (38) 
J 2tt L J 
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where the correction to the lepton self-energy is 

(39) 

In contrast to Eq. (120]) the propagator iS^^ in this expression is not the tree propaga- 
tor, but the propagator including a one-loop self-energy correction. Hence, we must also 
allow for off-diagonal terms and obtain for the wave-function corrections to the neutrino 
Wightman functions: 



~' l &N{j(k') — ^Ni^Nij^Nj iSm^Nij^Nj ^Ni^Nij^Nj + ^Ni^Nij^Nj 



k> 



(40a) 

— iS^f(k ) = iSxii'tf'NijiSjfj — iiS^i^^iS^ — iS^i^^-iS^ + iSj^^^iSj 



>Nj 



(40b) 



The different signs are a consequence of the Feynman rules on the CTP, and all functions 
on the right-hand side carry momentum argument k' as indicated after the bar. 

In order to reduce the number of terms that we need to write, we assume that 
M 2 Mi, such that the density of N 2 can be neglected at the time of leptogenesis. In 
that case, the off-diagonal components reduce to 

— i^ATUJ = 1 Sn1^N12 1 S'n2 ~ ^-^Nl^ N12^^N2 5 (41a) 

— 1^12 = ^N1^N12^N2 ~ ^N1^N12^N2 i (41b) 

~ i^21 = ~~ i^2^Ar21^ATl + ' 1 &N2^N2\ 1 S'n\ > (41c) 

~~ i^21 = ~' l Slv2 1 ^ N21 1 S'ni + ^N2^N21^N1 • (41d) 

As long as |M 2 — Mi| 3> 1^1,1^2, where IVi are the total decay rates of N iy we can 
straightforwardly extend the present analysis to situations where M 2 3> Mi does not 
hold by simply adding the corresponding contributions that arise from a non-vanishing 
density of iV 2 , i.e. the asymmetry produced in N 2 decays. The reason for this is that 
not both of the neutrinos in the wave-function diagram can be on-shell simultaneously, 
provided their mass terms are not degenerate. Because of this we do not further simplify 
iSj/2 to — i/M 2 at this point. 



4.2 KMS and Real Intermediate State Subtraction 

It follows from Eq. ( )38|) that vanishes provided both ^2J f<,> and Sf^ satisfy the 
KMS relation. The correction ^2™ f<,> in turn satisfies KMS provided it also holds for 
iStf i ^" > , see Eq. ( 1391) . It is instructive to show this property, 

iSjfcW = -ef"«iS$}(k) (42) 
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explicitly, because it illustrates how the RIS subtraction is realised within the CTP 
framework. 

To show the validity of the KMS relation for the wave-function correction ( H2|) ex- 
plicitly, we define 

K l3 =iS^{k) + e^iStUk)- (43) 

The claim is /Qj = when substituting equilibrium distributions. We show this on the 
example /C12, from which it will be clear how to generalize to all components of /Qj. 
Applying KMS to % N \ 2 and S^'i , we find 

/C12 = — iSj^i (^ m2 + ?7V12) l ^N2 + ^^Nl^N12^^N2 + ^N1^N12^N2 ■ (44) 

In addition, we have used that SjV2 is evaluated here only when N 2 is off-shell, such 
that we can set Sjf 2 = —Sjf 2 . To further simplify this expression, we make use of the 
following identities: 

S m + S m = S^i + , (45) 

% T m 3 + %% = + ■ (46) 

Note that the first identity, applied to the heavy neutrino N 2 , recovers the property 
Sjf 2 = —Sjf 2 , when evaluated for off-shell momenta. It follows that 

/C12 = -iS^ii $Ni2 + ^nu) '^m + K^m + ' 1 ^ni)^$ni2^N2 = > (47) 

where we have once more applied the KMS relation to obtain the last equality. An imme- 
diate consequence of this result is that no asymmetry is produced in thermal equilibrium, 
as required. 

Besides verifying the KMS relation for and consequently also for and ^ ' , 
the above calculation shows how the RIS subtraction emerges in the CTP framework. 
To see this, we neglect the finite density corrections in S £ ' implying that we set f\ = 
= /</> = /</> = in the propagators of the T and T type. Then, in order to see which 
process a certain contribution in the collision term corresponds to, we cut through the 
propagators of the <, > type, while leaving the propagators of the T, T type intact. 
We can then identify the terms in Eq. (jUJ) that involve ^jjj and with decays 

and inverse decays. In the Feynman diagram representation of that is given in 
Figure (TJ this corresponds to the cut indicated by the dashed double line, which gives 
the interference term between the tree-level decay amplitudes and the wave-function 
corrections to the decay amplitudes. The terms in Eq. (JUJ) that involve and Sjf^ 

correspond to scatterings. Since Sj^ and Sj fl are added together, precisely the pole 
contribution, when N± is on-shell, is isolated from the scatterings, which takes the effect 
of RIS subtraction. In the Feynman diagram in Figure [TJ this corresponds to the cut 
indicated by the solid double lines. It represents the the squared amplitudes of the 



13 




Figure 1: Diagrammatic representation of the lepton- number violating contribution 
to ^J f<,> . The solid lines with arrows represent the lepton £, the solid lines without 
arrows the neutrinos A^ 12 and the dashed lines with arrows the Higgs boson <p. The solid 
(light grey/orange) double line represents the cut that corresponds to the subtracted real 
intermediate states when finite density corrections to CP violation are neglected. The 
dashed (dark-grey/blue) double line represents the cut that corresponds to decays and 
inverse decays when finite density corrections are neglected. 

scatterings of leptons and Higgs bosons through N\ and N2 as well as the interference 
terms between these, which are crucial for the cancellation of CP- violating effects in 
equilibrium. Within the full CTP framework, and the T, T-type propagators also 
acquire contributions from real particles in the plasma, such that the clear separation in 
decay and scattering processes through the cuts in Figure [D no longer persists. 

In order to see the role of the violation of time reversal symmetry, we note that for 
the vanishing of /Qj, we make use of the KMS relation in the form of 

iS>#)iS<(fc')iA^") = iS^iSnk'yiA^k") (48) 

for momenta satisfying k = k'+k". This equation is valid not only in thermal equilibrium, 
but it is also another way of expressing that the process iVj — > £<p and its inverse £<p — > Ni 
occur at the same rate, and that the same holds true for TVj — > £<p* and £<fi* — > TVj. 
Hence, the above discussion does not only show that no CP asymmetry is generated in 
equilibrium, but more generally that no asymmetry is generated in situations that are 
symmetric with respect to time reversal, as it is of course expected by the CPT theorem. 

4.3 CP Source 

We now derive the CP-violating source within the lepton collision term, that arises due 
to deviations of N± from equilibrium and a violation of the KMS relations. First, we 
note that because of Eq. fl6]) we can identify 

iSS^i = iSStf! = iSSj fl = iSSj fl = i5S N1 , (49) 
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where iSSp^ ,T ' T denotes the difference between the non-equilibrium neutrino propagators 
and their equilibrium versions. Using this the CP-source collision term ( )38|) reads 



Cf(k) 



dk d A k' d A k" 
2tt (2tt) 4 (2tt) 4 



(2it) A 8\k — k' — k") 



(50) 



x 



Y*Y 2 i6S N1 i[$ 



J N12 



I l ^N2 



Y-iY^'iS^ 



^N2l ^ 



AT21 



iSS 



Nl 



J k' 



iA<(-k")iSf(k) 



+ 



Y*Y 2 i6S m i[$ 



J N12 



fj< N12 I 



A" 2 



^1 ^2*1^2 ( ?AT2: 



A^l 



i5S 



Nl 



J fc' 



xiAj(-fc")iSf(*) 



In this expression, the terms involving the dispersive parts ^^ T ' T of cancel. To 

see this, first note that i$^ T = — i^^ T . Then, since 5S , tv 1 (/c , k) = 5Sm(— ko, k) and 
since for equilibrium distributions iA^(A;Q, k.')iS^(k' , k') = iA^(— k / )i5^ > (— &q, k'), the 
cancellation occurs when performing the integrals. We can therefore make the replace- 
ments 



and write the wave-function correction to the lepton collision integral as 



d 3 k 



c 



wf / 



d A k d A k! d A k" 
(2vr) 4 (2tt) 4 (2tt) 4 



(27r) 4 5 4 (fc — kf — k") 



(51) 



(52) 



x tr 



-F^Y^Sati- {^ N i2 - *$NVl) l ^N2 + ^l^lS^g (?A>21 ~~ ?AT2l) 



ATI 



x (iA^-fc")i^<(fc) -iA>(-fc>S>(fc) 



The spectral self-energy is 



d 3 fc' 



d 3 fc" 



_(27r)V(A;-A;'-A;") 



(27r) 3 2|k'| 2 (27r) 3 2|k" 
x sign(A;o) [F^'Pr + x [1 - f? (k') 

In further simplifying C™ f , we obtain another factor of sign(fc ) from 



iA<(-k")iSm-^(-k")iSf(k) 



sign(fc )=-sign(fc ( / ) ') 



(53) 

/; q ( k ")] • 

(54) 



-P^Pr (27r) 2 5(k 2 )S(k" 2 )sign(k ) (1 - fP(k) + /; q (k")) 
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Using this result in Eq. (152]) . the tree- level propagators, and substituting after performing 
the Dirac trace the definition ( 133]) . we obtain our final result for the CP- violating wave- 
function correction to the lepton collision term: 

swi = I (0 c " (k) (55) 

(27r) 4 5 4 (fc' - k - k") 



(27r)32|k| (2vr)32 v /k' 2 + Af 1 2 (27r)32|k 

M X M 2 
k' 2 - Ml 



2lm[Y 1 2 Yf]2k^ N (k')5f N (k') (l - /f (k) + f;\k")) 
M X M 2 f d 3 k' E JVM (k')E^(kO 



Ml -Ml) (2vr)32 v / k 7J TM| 



Note that at zero temperature, we find 



T=0 l07T 



Then, by comparison with the tree-level lepton collision term (1201) . we reproduce the 
zero-temperature result for the decay asymmetry 

f rwg - r^gt 1 imiify; 2 ] m x m 2 
T=0 r^ w + r m ^ t 8tt ir^a - Mf ' 1 j 

where T denotes the decay rate of the process indicated in the subscript. On the other 
hand, if finite-density effects are kept for the initial and final states but not in the 
self-energy loop, then in the final expression of Eq. (I55D one substitutes 



E^(k')E^(k') > fc'^(k') ^ 

16tt 

5 Vertex Contribution to the CP Asymmetry 
5.1 Vertex Correction to Ci 

The vertex correction contribution to the CP asymmetry arises from the two-loop self- 
energy diagram for the lepton propagator (see Figure [2]) : 

ir e ab (k) = - cdYfY 2 J P*}SU-P)C {P^StiP + k + q)P R \ t C^ (59) 

x iS^i-q^iAfi-p - k)iA b ;(-q - k) . 

Here, the transposition t acts on the Dirac indices, and i, j and the CTP indices c,d = ± 
are summed over. The charge conjugation arises as a consequence of lepton number 
violation and the reversed flow of the internal lepton in the self-energy diagram. 
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We again consider the case where M 2 3> M\ in order to reduce the number of terms 
to be considered by one half. Furthermore, we neglect the diagonal (i = j) contributions 
to the self-energy since they do not contribute to the CP asymmetry. This can be 
easily verified along the lines of the calculation presented in this Section. With these 
simplifications, we obtain 



WW 



d 4 p d A q 



(60) 



x 



{y* 2 y 2 2 



(2vr) 4 (2tt) 4 

iS> Nl (-p)C [P L iSj(p + k + q)P R ] t C^ N2 {-q)iA<{-p - fc)iAj(-g - k) 
\S T m {-p)C [P L iS<(p + k + q)P R f C^Sjj 2 (-q)iA<(-p - k)iA<{-q - k) 

iSU-p)C [PiiSf(p + k + q)P R f ChS^i-q^i-p - k)\A<{-q - k) 
iS% 2 {-p)C [P L iSj(p + k + q)P R } t C^S> 1 (-q)iAl(-p - k)iA<{-q - k) 



Y 2 Y 2 2 



The expression for [^^(k) is obtained through interchange of the CTP indices T -h- T 
and < <H- >. 



5.2 KMS and Real Intermediate State Subtraction 

Before proceeding with the calculation of the CP source, we again show that this con- 
tribution leads to a vanishing asymmetry in equilibrium, which in particular means that 
the KMS relation, ilp^ [k) + e /3fc °i^ < (/c) = 0, must hold when equilibrium distributions 
are inserted. To see explicitly how this arises, we consider the contributions oc Y{ 2 Y 2 2 . 
First, we define 



(61) 



(27r) 4 (27r) 4 

and note that by virtue of the KMS relation 
J{k,p,q) = Y* 2 Y 2 2 [iSj ll (-p) + iSl 1 (-p)]c[P L iS<(p + k + q)P R ] t CHSj f2 (-q) (62) 
xiA<(-p-k)iA<(-q-k) 
- YfYzHS^i-p^ [P L iSj(p + k + q)P R Al(-q - k) 

+ P L iSf(p + k + q)P R Al(-q - k)J ' C^S T N2 {-q)\A<{-p - k) . 

Further simplification is achieved when deriving from the explicit expressions (jl]) and ([8]) 
that for sign(r ) = sign(so), 



i^ T (r)iAl( S ) + LS?(r)iAl(a) = iS<(r)iA<(s) + iS>(r)iA>(s) . 



(63) 
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Figure 2: Diagrammatic representation of ^J < ' > . The solid lines with arrows represent 
the lepton £, the solid lines without arrows the neutrinos N\ )2 and the dashed lines 
with arrows the Higgs boson <fi. The solid (light-grey/orange) double line represents 
the cut that corresponds to the subtracted real intermediate states when finite density 
corrections to CP violation are neglected. The dashed (dark-grey/blue) double line 
represents the cut that corresponds to decays and inverse decays when finite density 
corrections are neglected. 



We find 

J(k,p,q) = ^^[^(-pHi^-^ (64) 
xiA<(-p-k)iA<(-q-k) 

- Y* 2 Y 2 2 \S> ni (- P )c[p^S<( V + k + q)P R A<(-q - k) 

+P L iS>(p + k + g)P R A>(-g - k)] ' ChS T N2 {-q)\A<{- V - k) 
= 0, 

where for the last equality, we have again used the KMS relation. Likewise, one can 
show the vanishing of the terms oc Yi 2 Y 2 * 2 in equilibrium. 

Again, we can relate this verification of the KMS relation to the standard procedure 
of RIS subtraction. As in the wave- function contribution, at zero temperature the <, > 
type Green functions correspond to cut propagators. The terms in (16"U|) involving S^'i 
and thus correspond to the cut indicated by the dashed double line in Figure [2] 

that represents the interference between the tree-level and one-loop contribution to the 
N\ decay amplitude. The terms involving S N \ on the other hand correspond to N\ 
mediated scatterings, as indicated by the solid cut line in Figure El Note that here 
the cut corresponds to the interference between s-channel and t-channel diagrams, while 
the central cuts through the wave-function correction yield the squares of the s-channel 
and t-channel contributions to the scatterings. This complementarity is crucial for the 
separate cancellation of the asymmetries from the wave function and the vertex correction 
in equilibrium. 
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5.3 CP Source 

We now extract the CP source in the lepton collision term using the same simplifications 
as in the calculation of the wave-function correction. In particular, the deviation of N\ 
from equilibrium is parametrised again by SSm- We obtain 

c;(k) = / fr tr I^WW - WAQusrik)] (65) 

^(04 * *[{ ~ YfY 2 H5S N1 (-p)C[iSj(p+ k + g)iAj(-g - fc) 

-LS?(p + fc + g)iA<(-g - fc)]'cti^ 2 (-g)iA<(-j9 - fc) 
+Y 1 2 Y* 2 iSj I2 (-p)C[iSj (p + fc + g)iAj(-g - k) 

-iSfip + k + g)iA<(-g - *0]*C7tiA^(-g)iA<(-j, - k)]iS<{k) 

- Y^YfiSU-P^ [iSf(p + k + g)iAj(-g - k) 

-iS>(p + k + g)iA>(-g - k^ChSS^i-q^i-p - k) 
YfY 2 2 i5S Nl {-p)C[iSj{p + k + g)iAj(-g - fc) 
-LSf (p + fc + g)iA>(-g - A0]'C7ti5* (-g)iAj(-p - fc)) LSf (*)' 



where i<S , ^ 2 = — f° r off-shell neutrinos iV 2 has been used. Inspection of the Dirac 
structure reveals that Sjvi,2 only contribute through terms proportional to M 12 , which 
is familiar from the standard calculation of the decay asymmetry |T] . Finally, similar to 
the replacement (|5"T1) we can substitute here 

S T /A T f - S>A> +1 (S<A< - S>A>) , (66) 



SJ> T A T / - S<A< -+ -\ (S<A< - S>A>) 



under the integrals. 

The expression ( 165]) closely resembles the result for the wave-function contribution, 
except that iS^ 2 now depends on the loop momentum q. In the limit where M 2 S> Mi, 
the N 2 propagator can be expanded in powers of fc 2 /M| where k = (p, k), showing that 
in this limit the vertex correction is just one half of the self-energy contribution, where 
the factor one half arises since there is no summation over SU(2) indices within the loop. 

Inserting the explicit expressions for the lepton and Higgs boson Green functions, we 
obtain our final result for the vertex correction to the CP asymmetry. To this end we 
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introduce 

A « = J tM^^W^^ ~ * - ^ "> (P Zt- Ml ^ 
x [1 - /,(k') + U(k")\ . 

This has to be integrated over the initial and final state phase space in the collision term, 
such that we additionally define 

x [1 - /i(p') + Mp")1 ■ 
In terms of these quantities, the lepton collision term can be expressed as 

SV = j ^ C ^P') = 4Im [ y i 2 ^ 2 ] / ^Sf m (k)V(k,M 2 ). (69) 

5.4 Source term in the strongly hierarchical limit 

For the strongly hierarchical case, M 2 3> Mi, which we will focus on for the numerical 
analysis, the vertex contribution is simply one half of the self energy contribution. To 
see this, we note that in the limit M 2 3> Mi, 

MiM 2 Mi . . 
— 2 "> ~TT ( 7 °) 



(k f - p") 2 - Mf M 2 
in Eq. (I67p . which implies 



Mi E M (A) , „ . , _ M x E^(k')E^(k') 



r,(fc,/,M 2 ) = --^^^ and V (k,M 2 ) = - 1 -± ^ ' a (71) 
M 2 g w M 2 g z w 



Inserting this into (1691) and adding the wave-function contribution (155|) the source term 
for the lepton asymmetry in Eq. (13TI) is given by 

^M 2 »M! - - X 41m y x y 2 J -— / — 0/M k ) , 

2 V M 2 J J (27r) 3 2Vk /2 + Mf 9w 

(72) 

with Sjv^(k') given in Eq. (155]) . 



6 Numerical Estimates of Finite Density Effects 

In this Section we solve the kinetic equation for the final lepton asymmetry with CP- 
violating collision terms derived in Sections |3]- [5j We compare our results with all finite 
density effects included to the case where these effects are included in the initial and 
final state phase-space distributions, but not in the quantum corrections. 
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6.1 Effective CP-violating parameter 

To begin with, however, we present a rough estimate for the total magnitude of CP 
violation in the lepton collision term. This estimate is provided through an effective 
CP- violating parameter, which we define as: 

4-{n t -n t ) - W 
e cS = ~^r^ > (73) 

where n^m — J (f^w/4ivi(k) are the comoving number densities of lepton and neutrino. 
The tree- level (washout) contribution W, see Eq. ( )37l) . is subtracted from this definition, 
so that the numerator involves only the CP source, which is comprised of self-energy 
and vertex type corrections to the lepton collision term, given by Eqs. ( 1551 1691) . The 
CP- violating parameter thus also breaks into the corresponding two parts. For the 
self-energy part we get by using Eqs. (I55I) : 

£Wi = £Wi 16n Jd\k\^ m (k)^ Nlll (k)6f m (k) 
£ ° S £T ~° 9 W fd\k\^k^ N1 (k)5f N1 ^) ' 

where e^L is the CP- violating parameter in the vacuum, given by Eq. fl57|) . We see right 
away, that Eq. ( 1741) indeed reduces to e^L Q in the limit where the finite density effects in 
the loop are neglected in which case one factor of S^(fc) in the numerator is replaced by 
gwk^ / (IQtt) . We will not compute the vertex correction part here in full generality, but 
we note that in the strongly hierarchical case, where M2 S> Mi, the vertex CP- source 
reduces to one-half of the self energy CP-source, as shown in the previous Section. The 
same is then obviously true for the corresponding CP- violating parameter el s , so that 
the ratio £ e ff/£T=o is the same for the total and the wave-function contribution alone. 

To compute e™l requires the knowledge of the neutrino distribution function /jvi(k) 
for each instant of time (or temperature in the expanding Universe), which is obtained by 
solving the tree-level Boltzmann equation ( B1J) . Before doing this we make an estimate for 
ejg by factoring out the nontrivial time-dependence in Eq. (I74p . This can be achieved 
by assuming that the decaying heavy neutrino N\ remains in kinetic equilibrium and 
making an ansatz with a pseudo-chemical potential Uni- 

/jvi(k) = 1 . (75) 

exp[/3( v /k 2 + M 1 2 - m )] + l 

For small deviations from thermal equilibrium, this can be expanded in linear order of 
fi N1 as 

^i(k) = ^ q 1 (k)(l-/^ q 1 (k))^. (76) 

Using this ansatz the pseudo-chemical potential /xjvi involving the time dependence can- 
cels out in the ratio in Eq. (1711) . and we get 

£Wi = £W{ 167T Mkjjg s^aqsw*) /"i( k ) g - gW) m 

£cS £t=q 9 W j d\k\^ k^k) aw (i - am 
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To evaluate this expression, we see that the phase space integrals in E^ ri (A;), Eq. f )33|) . 
can be computed to give 

*<*> = S 7l @'r)' (78a) 

T i ,~ 9 W T 2 fko fko |k|\ M x 2 /*b |k|\\ fe* 

with 

2(2/0+2/) 

4(y ,y) = / ^^(i-^Y + ^irr) . (79) 
3(2/0-3/) 

which have analytical expressions 
f e yo+y — \\ 

^o,y)=ln^^ 7 - T J-y, (80) 
1 /l + e^\ 1 /l + e^\ 

A(j/o, y) =3(2/0 + y) in I -yn+« j - 2 ^° ~ ^ ln I i_ e "T y j ^ 

+ Li 2 ( -e 2 J + Li 2 f e 2 j _ Li 2 f-e 2 j - Li 2 f e 2 » 

where Li 2 denotes dilogarithm. The remaining one- dimensional momentum integrations 
in Eq. (177]) can then easily be computed numerically We plot the resulting ratio / e^L 
(equal to e c h/et=o in the strongly hierarchical limit) as a function of z = Mi/T in 
Figure [3j We find that the medium effects tend to enhance the effective asymmetry, 
in particular, the Bose enhancement of the Higgs particle dominates over the Fermi 
suppression of the lepton. 



6.2 Numerical solution of the Boltzmann equations 

We now scrutinise how the enhancement of the asymmetry is reflected in the solution 
of the Boltzmann equations ( )30|) . ( 13TT) . For our numerical analysis we consider the 
asymmetries in the limit of strong hierarchy: M 2 3> Mi, when the total CP source in 
Eq. ( 15T]) is given by Eq. ( 172]) . The vector can be evaluated according to Eqs. (178"]) 
when substituting for T the "comoving" temperature T com = a(rj)T. The washout term 
on the right-hand side of Eq. (131 j) has been given in Eq. ( 157]) . 

In the parametric region of strong washout, the neutrino density and also the lep- 
ton asymmetry will eventually become independent of the initial conditions, while this 
is not the case in the weak-washout domain. When dropping the assumption that the 
right-handed neutrinos are in equilibrium [161 HZj , we can integrate the equations for 
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Figure 3: The ratio of the CP-violating parameter when taking effects of finite density 
in the loop into account over the CP-violating parameter at T = 0. 

the comoving modes k com of the neutrino number distribution separately, which is the 
method that we pursue here. Thus, to solve the time-evolution equation for the lep- 
ton asymmetry, we first calculate the out-of-equilibrium neutrino distribution by solving 
Eq. (j30j) with the decay and inverse decay term P(k com ) given by Eq. fl34|) . This solu- 
tion is then inserted into the source term. The expression flM|) carries dependence on 
conformal time r\ through the dependence on k com0 = \/k-l om + a^rfMf, and the depen- 
dence of E^(/c com ) and /^\(k com ) on /c com o- Thus a numerical solution of Eq. ( 130]) for 
the out-of-equilibrium neutrino distribution is necessary. 

In Figure HI we show solutions to the Boltzmann equations in the intermediate regime 
between weak and strong washout. What is shown is the asymmetry summed over spin 
and SU(2) degree of freedoms and normalized to the entropy density, 

Y e = 2g w -^ = 2g w ' — '-, (82) 

S n I ^ 

45 i>* com 

where the last form is convenient since we calculate directly the comoving particle den- 
sities. We assume an effective number of relativistic degrees of freedom according to the 
Standard Model (g* = 106.75). In the Figure, we provide results both with (solid lines) 
and without (dashed lines) the additional finite density enhancement in the loops. We 
choose the parameters M l = 10 13 GeV, M 2 = 10 15 GeV, \Y 1 \=2x 10" 2 , \Y 2 \ = 10" 1 and 
a maximal CP phase such that ImfYj 2 !^* 2 ] is purely imaginary. Besides, we consider two 
different initial conditions, namely vanishing (light grey/red lines) and thermal (dark 
grey/blue lines) initial distributions for Aq. Recall also that the parameter z = Mi/T 
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Figure 4: The absolute value of lepton-to-entropy ratio Yi over z = M\jT . The param- 
eters are M x = 10 13 GeV, M 2 = 10 15 GeV, Y± = 2 x 10" 2 , Y 2 = lO^ 1 and a maximal CP 
phase. Both, thermal initial conditions (dark grey/blue) and vanishing initial conditions 
(light grey/red) for Ni are chosen. The solid lines correspond to solutions where finite 
density corrections in the loop are taken into account, the dotted lines to solutions where 
these are omitted. 

is proportional to rj. We observe from the Figure that there is a sizable effect of order 
(10 — 20)% from the finite density terms in the loops surviving at late times, since the 
asymmetry is partly produced at early times, when M\ is not much larger than T and 
quantum statistical corrections are of importance. This situation no longer persists at 
larger coupling, Y 1 = 5 x 10~ 2 , corresponding to strong washout. We display the result 
in the upper panel of Figure EJ Here the final asymmetry is determined at late times 
when Mi ^> T, such that quantum statistical factors have no sizable impact. 

Finally, we consider a scenario of weak washout (here we choose Y\ = 10 -2 ). The re- 
sulting asymmetries are displayed in the lower panel of Figure [5) In the case of vanishing 
initial conditions for N±, we find the impact of the finite-density effects in the loop in this 
scenario quite drastic: When they are taken into account, the asymmetry produced at 
early times when the iVVabundance is below its equilibrium value becomes larger, due to 
the Bose enhancement of the Higgs particles. This "wrong-sign" asymmetry is then too 
large to be cancelled by the "right-sign" asymmetry produced at later times, when the 
TVj-abundance is above its equilibrium value. Not only the temperature dependence of 
the washout, but also the temperature-dependence of the effective CP-violating param- 
eter is crucial for determining the final lepton asymmetry. This feature does not occur 
when the loop enhancement is neglected, since in that case, the effective CP asymmetry 
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Figure 5: The absolute value of lepton-to-entropy ratio over z = M\jT . The param- 
eters are M x = 10 13 GeV, M 2 = 10 15 GeV, Y 1 = 5 x 1(T 2 , F 2 = KT 1 and a maximal CP 
phase in the upper panel (strong washout) and Yi = 10~ 2 , Y 2 = 10 _1 in the lower panel 
(weak washout). Both, thermal initial conditions (dark grey/blue) and vanishing initial 
conditions (light grey/red) for Ni are chosen. The solid lines correspond to solutions 
where finite density corrections in the loop are taken into account, the dotted lines to 
solutions where these are omitted. 
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is temperature independent and the resulting asymmetry is always of the "right- sign". 
Therefore, neglecting the finite density effects in the loop in weak washout scenarios is 
in general not justifiable and may under certain circumstances even lead to an error in 
the predicted sign of the asymmetry. 

7 Conclusions 

We presented a full derivation and numerical solution of the kinetic equations describing 
leptogenesis in the CTP framework. The results extend existing approaches since all 
quantum statistical factors are accounted for systematically for the first time in a realistic 
model of leptogenesis. We have put particular emphasis on establishing a connection to 
the conventional Boltzmann approach and elucidated how RIS subtraction is realised in 
the CTP formalism. The numerical analysis shows that finite density corrections do not 
play a significant role in the strong-washout regime while they may have a sizable impact 
in weak-washout scenarios. 

The detailed analysis of the simple model of leptogenesis presented in this work 
should be seen as a first step toward applications of the CTP formalism to situations 
where the Boltzmann framework needs to be modified. One may think of flavour ef- 
fects [231 12H 1251 [26] , the inclusion of additional thermal effects [15] or resonant leptogen- 
esis J27J El ESI [21] in the limit of degenerate right-handed neutrino masses. Examples of 
approaches to these more advanced problems using the CTP formalism may already be 
found the literature [TOj [Til [29] . However, so far there has been no complete derivation 
and solution of the equations describing the standard scenario of leptogenesis within the 
CTP framework. The present work may be viewed as a self-consistent leading-order ap- 
proximation in this framework, which serves as a starting point for investigating further 
corrections and variations. 

Note added 

While this paper was in writing, Ref. [30J appeared which also studies modifications of 
an effective CP- violating parameter in a realistic leptogenesis model due to finite-density 
effects in loops. 
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